KEK-TH-1021 
hep-th/0507029 



in 
o 
o 

(N 



X 



Thermodynamic Behavior of Fuzzy Membranes 
in PP-Wave Matrix Model 



Hyeonjoon Shin^ and Kentaroh Yoshida 



*b 



t BK 21 Physics Research Division and Institute of Basic Science 

Sungkyunkwan University, Suwon 440-746, South Korea 
> • 

Q\ ' * Theory Division, High Energy Accelerator Research Organization (KEK), 

(N ! 

O ■ Tsukuba, Ibaraki 305-0801, Japan. 

X> 

o 
in 
o 



a hshin@newton.skku.ac.kr 6 kyoshida@post. kek.jp 



Abstract 

We discuss a two-body interaction of membrane fuzzy spheres in a pp-wave matrix model 
at finite temperature by considering a fuzzy sphere rotates with a constant radius r around 
the other one sitting at the origin in the SO (6) symmetric space. This system of two fuzzy 
spheres is supersymmetric at zero temperature and there is no interaction between them. Once 
the system is coupled to the heat bath, supersymmetries are completely broken and non-trivial 
interaction appears. We numerically show that the potential between fuzzy spheres is attractive 
and so the rotating fuzzy sphere tends to fall into the origin. The analytic formula of the free 
energy is also evaluated in the large iV limit. It is well approximated by a polylog-function. 
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1 Introduction and Summary 



The basic degrees of freedom of string theory and M-theory are fully encoded in matrix models 
[1-3] . The matrix models are believed to give non-perturbative formulations of string theory and 
M-theory. In particular, the BFSS matrix model is supersymmetric matrix quantum mechanics. 
It is believed to describe a discrete light-cone quantization of M-theory. It also describes the 
low-energy dynamics of iV DO-branes of type IIA superstring theory [4]. In addition it is a 
matrix regularization of the light-cone action for the supermembrane in eleven dimensions [5]. 

Matrix model thermodynamics is also interesting in relation to black hole physics. One 
motivation for understanding the behavior of the BFSS matrix model at finite temperature 
comes from the conjecture that their finite temperature states are related to black hole states 
of type IIA supergravity [6,7]. This idea was studied in a series of papers by Kabat, Lifschytz 
and Lowe [8]. The BFSS matrix model at finite temperature is also investigated in [9]. Some 
features are discussed in [10] (For a review of brane thermodynamics, see [11]). 

By the way, a matrix model on a pp-wave background was proposed by Berenstein-Maldacena- 
Nastase (BMN) [12], and it has been intensively studied. The background of this matrix model 
is given by the maximally supersymmetric pp-wave background [13]: 

ds 2 = -2dz + cfe-- ^T(|)V) 2 + E^ (1-1) 

\t=l a=4 / 7=1 

-F+123 = A* • 

The action of the matrix model on this background .Spp consists of two parts as follows:* 

(1-2) 

^DtX'DtX 1 + ^([X I ,X J )) 2 + i Q^D t Q-RQ^ I [e,X I ] , (1.3) 

where the indices of the transverse nine-dimensional space are /, J = 1, . . . , 9 and R is the 
radius of the circle compactified along x~ . All degrees of freedom are N x N Hermitian 
matrices and the covariant derivative D t with the gauge field A is defined by D t = d t — i[A, } . 
This matrix model is closely related to the supermembrane theory on the pp-wave background 
via the matrix regularization [5] (For works in this direction see [14,15]). This matrix model 
* Hereafter we will rescale the gauge field and parameters as A — > RA , t — > t/R, fi — ► R/j, . 
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has a supersymmetric fuzzy sphere solution (which is called "giant graviton" ) due to the Myers 
effects [16], because the constant 4- form flux is equipped with. This fuzzy sphere solution is 
given by 

-^sphere = 3^ ' (1-4) 

where J 1 satisfies the SU(2) algebra [J 1 , J J ] = ie^ k J k . The classical energy of this solution is 
zero and hence the fuzzy sphere can appear in classical vacua without loss of energy. Namely, 
the classical vacua of the pp-wave matrix model are enriched with fuzzy spheres, and it may be 
interesting to look deeper into the dynamics of fuzzy sphere solution. Then the fuzzy sphere 
solution ^Q phere preserves the full 16 dynamical supersymmetries of the pp-wave and hence is 
1/2-BPS object. We note that actually there is another fuzzy sphere solution of the form 
but such solution does not have quantum stability and is thus non-BPS object [17] (For other 
classical solutions, see [18-21]). In [22] the interaction potential of membrane fuzzy spheres 
was shown to be the l/r 7 -type. 

In this letter we discuss a thermal interaction between fuzzy membrane solutions in the 
pp-wave matrix model. The two-body interaction of fuzzy spheres at zero temperature was 
studied in our previous work [23] (containing the extension of the work [24]) by considering the 
system that a fuzzy sphere rotates with a constant radius r around the other one sitting at 
the origin in the transverse six-dimensional space (Fig.^). In the case of zero temperature the 
interaction potential is zero basically because of the existence of the remaining supersymme- 
tries. In the finite temperature case, however, non-trivial interaction should appear since the 
supersymmetries are completely broken due to the thermal effect. Here, we will be interested 
in the evaluation of the thermal potential at finite temperature by heating the system depicted 
in Fig.QJ Firstly we consider the exact one-loop free energy by using the /i — > oo limit. In 
this computation we use the spectrum around the two-body system which has been obtained 
in [23] . By numerically plotting the free energy, we can see that the rotating fuzzy sphere tends 
to approach to the static fuzzy sphere. The attractive potential becomes stronger and stronger 
as temperature grows and the distance decreases. Furthermore, as the size of the matrix be- 
comes larger and larger, the potential becomes deeper and deeper. In addition we evaluate the 
analytical expression of the free energy by taking the large N limit. Then the discrete sum 
for the index of spin may be replaced by the integral for it. When we utilize the long-distance 
expansion, we find that the leading part of the free energy can be represented by a polylog 
function. The approximation of the free energy qualitatively well reproduces the behavior of 
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Fig. 1: The two-body system of fuzzy spheres. A fuzzy sphere rotates around the other one sitting at the 
origin in the transverse six-dimensional space. 



the exact free energy. 

This paper is organized as follows: In section2 we introduce the formula for the exact free 
energy of a pp-wave matrix model at finite temperature. In section 3 the free energy for the 
interaction between two fuzzy spheres is evaluated both in numerical and analytical methods. 

2 One-Loop Exact Free Energy at Finite Temperature 

We will briefly review the calculus of the one-loop free energy in the finite temperature case. 
By using the background field method and the /x — > oo limit, it can be exactly computed. 
We first decompose the variables as 

x I = s I + y I , e = o + #, (2.i) 

where B 1 are the classical background fields while Y 1 and \I/ are the quantum fluctuations 
around them. In some cases we need to consider the background of the gauge field (gauge field 
moduli) due to the finite temperature effect. In the case we consider here, however, we do not 
need to take account of it, as we will explain later. 

In order to perform the path integration, we take the background field gauge which is usually 
chosen in the matrix model calculation as 

D*A* l = D t A + i[B I 7 X I ] = 0. (2.2) 

Then the corresponding gauge-fixing Sqf and Faddeev-Popov ghost Sfp terms are given by 

S GF + S FP = JdtTr^-^D^A^-CdtDtC+lB 1 ,^ 1 ,^ . (2.3) 



Now by inserting the decomposition of the matrix fields (j2.1j) into the matrix model action, we 
get the gauge fixed plane-wave action S (= S pp + Sgf + Sfp) expanded around the background. 
The resulting acting is read as S = So + S2 + S 3 + S4 , where S n represents the action of order 
n with respect to the quantum fluctuations and, for each n, its expression is 



S 

S-2 



dtTi 
dtTi 



i(F 7 ) 2 - 2iB T [A, Y 1 ] + ^([B 1 , Y J ]f + [B 1 , B J ][Y ] , F J 



ifie ijh B l Y j Y k 



5', 



dtTr 



5, 



dtTr 



- ii 2 - A]) 2 + CC+ [B 1 , C\[B', C] 



iY T [A, Y 1 } - [A, B T ][A, Y 1 ] + [B 1 , Y J ][Y T , Y J ] + &[A, 

*V[*, v 1 ] - i|e^ fc F 4 r j r fe - c] + [b 1 , c][Y T , c] 



(2.4) 



For the justification of one-loop computation or the semi-classical analysis, it should be 
made clear that S3 and S 4 of Eq. (|2.4j) can be regarded as perturbations. For this purpose, 
following [14], we rescale the fluctuations and parameters as 

A^yT^A, Y'^y-^Y 1 , C^yT^C, C - yT^C , t - y~H . (2.5) 

Under this rescaling, the action S in the fuzzy sphere background becomes 

S = S 2 + /i~ 3/2 S 3 + y~ 3 S 4 , (2.6) 

where S2, S3 and S4 do not have y dependence. Now it is obvious that, in the large y limit, S3 
and S 4 can be treated as perturbations and the one- loop computation gives the sensible result. 
Note that the analysis in the S2 part is exact in the y — > 00 limit. We can calculate the exact 
spectra around an N- dimensional irreducible fuzzy sphere in the y — > 00 limit, by following the 
method in the work [14] (For the detail calculation, see [14,23]). 

Now we are interested in the finite temperature case where the system couples to the thermal 
bath. In order to consider the thermal system with temperature T, let us move to the Euclidean 
formulation via the Wick rotation t — > it, and compactify the Euclidean time direction with a 
periodicity (3 = 1/T. Note that T is a dimensionless parameter now because of the scaling of 



time variable t — > R l t. This compactification leads us to encounter the summation instead of 
the momentum integral. We can easily compute this summation by using the formulae 

In (nV + M' 2 ) = 2 In sinh M' (for bosons) , (2.7) 

n=intcgcr 

V In (nV + M' 2 ) = 2 In cosh M' (for fermions) , (2.8) 

n=half integer 

and the fact that the fuzzy sphere configuration under our consideration is supersymmetric. 
Here we should note that the ghost fields obey the periodic condition rather than anti-periodic 
condition. Hence its contribution gives the sinh but the integration of the ghost field gives the 
inverse sign in comparison to bosons. This fact ensures the cancellation among unphysical and 
ghost degrees of freedom. Thus the free energy F = —T In Z is represented by 

F = T N i ln (1 - e_/3A/l ) " ~ TN * ln (! + e^*) " 2TN c In (l - e~ pMc ) . (2.9) 

i£Y,A 

The symbols Y and A denote the bosonic fluctuations of X and A, respectively. The ^ and 
C denote the fermionic fluctuations of and C . The free energy usually contains the zero 
temperature part, but this part does not appear in the present case since the fuzzy sphere 
background is supersymmetric at zero temperature. 



3 Free Energy of Fuzzy Membrane Interaction 

We are interested in the configuration of classical solution in which a fuzzy sphere rotates with 
a constant radius r around the other fuzzy sphere (see Fig. [I}. This system is described by 

. 4) = f4) (< = 1,2,3; 8=1,2), 
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B m = 


r cos | 


v6 / 



Bfi) = Tsin (J^tj ljViXiVi , otherwise = . (3.1) 



L NlxNl , - K±J ■ 

The B 1 , x are iV s x N s matrices. We take B 1 as N x N matrices and then N = Ni + N 2 . Let 
us concentrate on the N\ ^ N 2 case so that we do not have to take account of the gauge field 
moduli considered in [26-28] . If two fuzzy spheres of the same size are considered then we may 
have the gauge field moduli. This gauge field moduli plays an important role in the study of 
the Hagedorn transition in the transverse M5-brane vacua [29]. 

This system (J3.1j) was proposed in our previous work [23] and it was shown to be supersym- 
metric. That is, the interaction potential between them vanishes (one-loop flatness) [23] as well 
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Tab. 1: The spectr urn of the interaction part between two fuzzy spheres. 



as the quantum fluctuations around each of the fuzzy spheres [14,17,23]. This result comes from 
the supersymmetries of the theory. However, the supersymmetries are completely broken down 
and so non-trivial potential should appear. In our previous work [25] we studied the stability of 
the fuzzy sphere > in the case of finite temperature. In particular, we compared the free energies 
between the fuzzy sphere vacuum (X 1 = | J 1 (i = 1, 2, 3)) and the trivial one (X 1 = 0) , where 
all the fuzzy spheres are considered to be located at the origin in the transverse six- dimensional 
space. Then we are now interested in the interaction potential between the two fuzzy spheres 
being at a distance in the transverse space and will consider the off-diagonal fluctuations which 
are responsible for the interaction. 

In the computation of the partition function for off-diagonal fluctuations, we should diag- 
t Stability of fuzzy spheres is discussed in the IIB matrix model with Chern-Simons term [30] and mass 
term [31]. 



6 



onalize the quadratic action. After the diagonalization, we see the spectrum. The resulting 
spectrum is summarized in Tab.[T] By using the spectrum the free energy is expressed as 

|(iVi+Ar 2 )-2 i(ATi+JV 2 ) 

\(3F = (2j + l)ln(l-eW^^ 2 )+ £ (2j + 1) ln(l - eW^' 2 ) 

j=||JVi-JV 2 | j=||JVi-JV a |+l 
|(JVi+JV 2 )-l 

+ E 4(2j + 1) ln(l - eW^+S^^) 

J=i|JV!-JV 2 | 
i (AT 1+A r 2 )_i 

+ E + 

j=^\N 1 -N 2 \ 
I(JVi+JV 2 )-f 

E 2 ( 2 ^' + x ) 

J=i|iVi-iV 2 |-i 
|(ATi+iV 2 )-i 

E 2 ( 2 ^' + x ) 

j=i|JVi-JV 2 |+i 

Here the unphysical degrees of freedom are canceled out. That is, the above expression has no 
degrees of freedom for fuzzy sphere rotation, gauge field and ghosts. 

We will consider the free energy of the fuzzy membrane interaction in both numerical and 
analytical methods below. To begin with, we numerically study the behavior of the free energy. 
Then we evaluate the analytical expression of the free energy by using the large N limit. 

3. 1 Numerical Study of the Free Energy 

The expression of the free energy for the fuzzy membrane interaction is quite complicated and 
so it is difficult to study analytically the behavior of the free energy without any approximation. 
We will analyze the analytic behavior of the free energy in the next subsection by using the 
large N limit and the long-range expansion. Here let us however investigate numerically the 
behavior of the free energy. In Fig.El we give the numerical plot of the free energy for some 
cases. From the results we can see that the interaction between two fuzzy spheres is attractive. 

3.2 Analytic Behavior of the Free Energy in Large N Limit 

In the previous subsection we presented the numerical study of the thermodynamic behavior 
of the free energy for the interaction between two fuzzy spheres. From now on we will evaluate 



ln(l - eW^' 2 ) + ln(l - e W^^> 2 ) 



ln(l + e~^ r2+ ^ +1)2 ) + ln(l + e ~^M^ 



In 



(l + e-V^^+D 2 ) 



+ ln(l 



(3.2) 
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i) Ni = 3 and N 2 = 2 ii) iVi = 20 and 7V 2 = 10 

Fig. 2: The numerical plots of the free energy F for the two-body interaction. The F is plotted in the vertical 
axis with respect to the distance r and the temperature T . The left one is for the Ni = 3 and N 2 = 2 , and the 
right is for Ni = 20 and 7V 2 = 10 . 



the analytical expression by considering the large N limit. In such situations we can replace 
the discrete summation for the spin j by an integral with a continuum variable and so it is 
possible to evaluate an analytical expression for the free energy. For example, let us focus upon 
a part of the SO (3) bosons, <x, m . We consider the large N limit: 

JVi, N 2 > 1 with N x -N 2 = c, (3.3) 

where c is a non- vanishing constant. Then the summation for aj m may be replaced as follows: 

i(Ar 1+ Ar 2 )-2 N 



(2j + 1) In (l - e~^ r2+ ^ [l+1)2 \ — / dx (2x + 1) ln(l - e" 



(3.4) 



j=±\N 1 -N 2 \ 



In addition, in order to evaluate the above integral analytically, we utilize the long-distance 
expansion (large r limit) and expand the square root as 



r2 + l(x + l)^r + T i-(x + l) 2 + (9(l) 



By using the Taylor expansion, 



V + \ y2 + l y3 + " ' 
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o -i 

5 10 15 20 



Fig. 3: The error function. 



the free energy for a jm can be evaluated as 

oo 

Fa = 36T V^^|e-^(^) (18r2+1+27V+7v2) -e-^(^) (72r2+4+4|c|+|c|2) } 



n=l 



-3V2^£ erf '-^J^ -erf 



n=l 

oo 



iV+ 1 /2n/3\ , / \c\ + 2 /2n/3\ e 



6 V r \ 12 V t ) \ n 3 / 2 



-nf3r 



-36TV£ie-* + 3V2^f:| l-e,f(^^ (3.7) 



71=1 71=1 



where we have used that N and r is sufficiently large so that erf(x) = 1 (x ^> 1) as shown 
in Fig-El but we have supposed that r ^> N . When we consider the low temperature region 
T 1 jr , the second summation in the expression (|3.7|) can be ignored. Then, using a polylog 
function, 

°° z n 

polylog(a, z) = ^ — . 

n=l 

we can write down the analytical expression of the free energy as follows: 

F = -36TV ■ polylog(2, e _/3r ) . (3.8) 

For the other bosonic sector we can evaluate the free energy in the similar way and the same 
result as ()3.8|) . Then the fermionic sector and the ghost sector the contributions are evaluated 

as 

F 9* -36T 2 r ■ polylog(2, -e _/Jr ) . (3.9) 
Thus the total free energy is given by 



-36 x 8T 2 r 



polylog(2, e"^) + polylog(2, -e~^ 



144TV • polylog(2, e" 2/3r ) , (3.10) 




r 



i) the approximated free energy in large N ii) the exact free energy for Ni = 50 and N2 = 30 

Fig. 4: A comparison of the approximated free energy in large N and the exact free energy In i) the 
approximated free energy is numerically plotted. In ii) the exact free energy in the same region as i) is shown 
in order to see the validity of the approximation. 

in the region: 

T < - , iV < r , r, N : sufficiently large . (3.11) 

r 

Here we have used the formula: 

polylog(2, z) + polylog(2, -z) = ^polylog(2, z 2 ) . 

In particular, polylog(2, 1) = 7r 2 /6 and polylog(2, 0) = . The numerical plot of ()3.10|) is shown 
in Fig.0] We see that the behavior of ()3.10|) is quantitatively similar to that of the exact free 
energy with large N values. In the low temperature region the large N free energy for the 
membrane fuzzy sphere interaction is well described by a polylog function. 

Note that we cannot naively evaluate the high-temperature behavior from ()3.7j) since the 
n = 00 contributes and so we cannot take the error function as erf(0)=0. 
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